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Abstract 

The statistical mechanics of ID and 2D Ginzburg-Landau systems is 
evaluated analytically, via the transfer matrix method, using an expres- 
sion of the ground state energy of the quartic anharmonic oscillator in 
an external field. In the 2D case, the critical temperature of the or- 
der/disorder phase transition is expressed as a Lambert function of the 
inverse inter-chain coupling constant. 

1 Introduction 

The physics of ID Ginzburg-Landau systems, described by a polynomial Hamil- 
tonian, with quartic anharmonicities, is of utmost interest, for several domains, 
from quantum field theory to structural phase transitions. 2- or 3D aggregates 
of such chains describe anisotropic anharmonic solids. Examples of this kind 
of 2D systems are the ultra-thin ferroelectric films, obtained quite recently pQ. 
In ultra-thin polymeric ferroelectric films, the strong dipols constituting the 
polymer chain allow the manipulation of the individual monomer by a local 
field, making this class of materials very attractive for molecular electronics 
[2][Cai2004]. 

Theoretically, one of the most popular approaches to the study of Ginzburg- 
Landau anisotropic systems is the transfer matrix method. With this method, 
the evaluation of the free energy of a (N + 1) - dimensional classical system is 
equivalent to the evaluation of the ground state energy of a quantum N - dimen- 
sional system, described by an neffective Schrodinger equation. The classical 
(TV + 1) - dimensional system and its quantum N - dimensional counterpart are 
sometimes termed dual systems [5] . Any progress in understanding the physics 
of such a system is mirrorede by a progress in understanding its dual system. 
In traditional approaches [J], [5], the solution of the quantum problem is used 
in order to find its classical counterpart. However, in some Monte Carlo simu- 
lations [6], the classical 2D system was studied in order to solve the quantum 
ID system. 
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Technically speaking, the transfer matrix approach, as developed by Scalapino 
and his co-workers [4] , [7] is mainly a convenient method of evaluating the par- 
tition sum, expressed as a functional integral. The first decade of investigations 
of this subject - until the early '80s - was dominated by significant analytical 
results. For the next two decades, the efforts were mainly focused on numeri- 
cal investigations. However, in the recent years, some exact results, or at least 
analytical approximations, obtained in the study of the quartic oscillator, show 
a change of this tendency. They allow us to give analytical responses to some 
important problems, which have been received, previously, only numerical solu- 
tions. 

To this trend belongs also the present paper. We use an analytical approx- 
imation of the ground state energy of a quartic oscillator in an external field, 
obtained by Van der Straeten and Naudts [8] , in order to find the exact statisti- 
cal mechanics of a chain of classical anharmonic oscillators, and the mean-field 
approximation of the statistical mechanics of a planar array of coupled chains. 

In fact, the formula proposed in [5] for the ground state energy of a quartic 
oscillator in an electric field is merely an ansatz than a rigorously obtained re- 
sult. Also, it depends only on the external field - any dependence of temperature 
of the effective mass of the oscillator, essential for the description of the statis- 
tical mechanics of the linear chain, is lost, being replaced by some numerical 
parameters. In our approach, both these drawbacks of the Van der Straeten- 
Naudts solution are removed. The temperature dependence is restored using 
a comparison with similar results obtained for the ID Ising model in external 
field. The adequacy of the ansatz is checked and confirmed by the physical 
character of its predictions and, more than this, by the fact that it gives results 
obtained using alternative methods. This is an interesting example of the fact 
that, while studying dual systems, a progress made in understanding one of 
them is mirrored in the understanding of the other. 

The outline of this paper is as follows. In section 2, we shall sketch the 
general frame of the problem, at least for sake of clearly defining the notations. 
In Section 3, we shall expose in some details the Van der Straeten-Naudts solu- 
tion, in order to properly understand its advantages and drowbacks. In Section 
4, this solution, with restored temperature depedence, is used in order to find 
the electrical polarization of the anharmonic chain. Section 5 is devoted to the 
physics of 2- and 3D arrays of such chains, treated in the mean-field approxima- 
tion. A new analytic formula, for the critical temperature of the order/disorder 
transition, as a function of the inter-chain coupling, is obtained. Its small cou- 
pling limit coincides with a formula obtained in [10 , confirmin the adequacy of 
the ansatz contained in the Van der Straeten-Naudts solution. The last section 
is devoted to conclusions. 
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2 Chains of anharmonic oscillators 



Basically, the physical system to be studied in this paper is a chain of classical 
anharmonic oscillators, described by the Hamiltonian: 
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Each oscillator is an atom; i,j indicate lattice sites (equilibrum positions of 
atoms) ; u% - displacements of displacing atoms with respect to some heavy ions 
or reference lattice. Periodic boundary conditions are assumed. The coefficient 
A is defined by the attractive interactions of the mobile atom with the reference 
lattice, B - by short-range repulsive interactions, and - by elastic interactions 
between displacing atoms. If (1) describes a lattice which is unstable against 
a displacive transition, A < 0, B > 0, Cij > 0. In this situation, the potential 
energy on site appearing in (1) has two minima, at 
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We shall replace the discrete Hamiltonian (1) with a continuum representa- 
tion: 
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jl locates an 



where I is the lattice spacing, Co - the sound velocity, and Xj 
atom in the continuum representation. 

The equilibrum thermodynamics of this ID model can be obtained from the 
classical partition function: 



Z = j VuDpe- m{u ' p) 



where, as usual, 



and 



Z p = (2vrfc B T) 



N/2 



(4) 



(5) 



Z u = e- N ^ 

with e - the ground state energy of the Schrodinger equation: 
1 



2m* du 2 



(6) 



(7) 
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The free energy F of the chain is: 



F = -k B T In Z. p Z u = F p + F u (9) 

but, in fact, only the "interaction term", F Ul produces interesting physical ef- 
fects. The free energy per particle is: 

fu = eo (10) 
In the presence of an electric field, a term 

-Equi = -pu % 

is added to the on-site energy (the notation used in the r.h.s. intends to com- 
promize with that used in [8]), and the quantum Hamiltonian (7) becomes: 

1 d2 A 2 B 4 
H„ = ;r + — u H u 4 - pu (11) 

9 2m* du 2 2 4 F V ' 

So, the average value of the displacement can be written as: 

W - M 

where eo (p) is the ground state energy of H q . Also, using any of the Hamilto- 
nians (7), (11), we can compute the average value of other powers of the order 
parameter: 

\ U) d\A\' \ / dB 



The physics of such a system is quite simple; we shall discuss it starting 
with the discrete variant, described by the Hamiltonian (1). In the absence of an 
external field, each atom oscillates near one of the two minima, given by (2); the 
average value of the displacement is zero, at any temperature. For a planar array 
of such chains - so, for a 2D system - the situation is quite different. The inter- 
chain interaction limits the effect of fluctuations, and the average displacement 
along a chain, (u), can have a non-zero value, if the temperature is lower than 
a critical value, T c . So, the sistem displays an order /disorder phase transition, 
similar to a magnetic/paramagnetic transition, in a magnetic system. Referring 
again to the ID system, an external electric field can produce a state with (u) ^ 
0. Also, the effect of inter-chain interactions, in a 2D system, can be treated, 
in the mean-field approximation, as an effective external field, able to trigger a 
phase transition. For the description of such a phase transition, it is important 
to evaluate the critical temperature as a function of microscopic quantities, 
mainly of inter-chain coupling. Also, quantities like the order parameter (u) , 
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the polarization and the electric susceptibility play a similar role to magnetic 
moment, magnetization and magnetic susceptibility. In fact, there is an evident 
connection between the atomic chain with deep on-site wells and the ID Ising 
model, at least for small external fields. 

In principle, the formulae given in this section, mainly (10), (12), can provide 
a quite complete description of the " exact" statistical mechanics of the ID sys- 
tem (if the eigenvalue eo is known " exactly" ) , and a mean-field variant of such 
a description, for the 2D system. However, as it is wellknown, the eigenvalues 
of (11) cannot be evaluated exactly, even for p = 0. Recently, an approximate 
analytical expression for eo (p) has been obtained in [3] . We shall take advantage 
of this formula in order to find the polarization, the electrical susceptibility and 
other properties of the chain of classical anharmonic oscillators, and of the 2D 
aggregates of such chains. The outputs of our calculations will be also a test of 
the adequacy of the expression of eo (p) given by [8] . 

In the next section, we shall briefly expose how the expression of eo (p) has 
been obtained, in order to properly understand its validity and its limitations. 

3 The quantum double well anharmonic oscilla- 
tor in an external field 

We shall briefly describe the results obtained recently by Van der Straeten and 
Naudts [8] on this subject. The Hamiltonian: 

p2 i 

H=—+Y J \Q i (13) 

is written in terms of annihilation and creation operators a, a + of the ordinary 
harmonic oscillator with mass m and frequency ujq. Defining 



r 



2 



mujQ 



(14) 



and using a formula for normal ordering derived in |11) . the Hamiltonian (1) 
takes the form: 



1=0 \ v ^/ k=0 j=0 J J I 

The expectation value of (15) in a certain state \ipt >, depending on a pa- 
rameter t, will be minimized 

A ^ t \ H \^)=0 (16) 
in order to obtain an equation which can determine the parameter t. 
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For a quartic oscillator in an external field, Ai = — p, X2 = A3 = 0, A4 = 
j, 1 = 4:, and (16) becomes: 

(2t + 1) (2X2^ ~ ^) + 6A 4 r^ (2t 2 + 2t + l) = (17) 

If f is known - or, equivalently, if the choice of the state vector \ipt > has 
been done - eq. (17) will give the value of tq; this cubic equation has indeed one 
real root, if the potential has two wells (a < 0) . The choice of t has no clear 
physical significance. The authors adopt the variant t = N/2, where N is the 
value at which the \ip t > basis is truncated, for reasons of rapid convergence 
of numerical calculations, to is also fixed numerically. The choice of a shallow 
well, with a = 2X 2 = —2, (3 = 4A 4 = 1, transforms (17) in an equation with 
numerical coefficients, which gives a numerical value for r . 

Finally, the ground state eigenenergy in the presence of an electric field, 
eo (p) , is, for small and moderate values of the field: 



eo {p) — eo (0) — |a|ptanhwp (18) 

where the coefficients a, oj are determined numerically. Their "exact" value 
depends on the truncation of basis functions \ipt > . In fact, this formula of eo (p) 
is not rigorously deduced, but is merely an interpolation between the "small" 
and "moderate" field cases. The corectness of this ansatz could be evaluated 
according to its success of giving reasonable predictions for the statistical physics 
of the Ginzburg-Landau systems, through the transfer matrix method. 
For large fields, 

e (p) ~A + B oP 4 / 3 (19) 

with Aq, Bo — numerical constants. A similar - but more general - result has 
been obtained by Bronzan and Sugar [12], for the energy levels of a potential 
V(x) = A\x + x 4 . It is clear that, for large fields (with Bronzan and Sugar's 
notations, for A\ > 15), the effect of the quadratic (harmonic) term of the 
potential is irrelevant, at least for the leading terms. 

Let us define what we mean here by a "large" field. In the absence of an 
external field, the particle oscillates in a symmetric, double-well potential: 

Vs (x) = ^x 2 + ^x 4 (20) 
In the presence of an electric field, one gets a total, asymmetric, potential: 

V as (x) = ^x 2 + ^x 4 - px (21) 

The "electric" term favors one of the two wells; this one becomes deeper, 
while the other one - more shallow, disappearing when the electric field is larger 
than a certain critical value p c - It is easy to find p c by analyzing the extremum 
points of V (l s(x), which are the roots of the equation: 
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dV as (x) 3 

: = px + ax — p = U (22) 

ax 

where a < 0, /3 > 0, p > 0. This equation may have one or three real roots, 
according to the sign of the discrininant, proportional to the expression: 

I |3 4 
M , P_ 

27/3 4 

which changes its sign at the critical value: 

/4H 3 \ 1/4 

If we presume a Landau-type temperature dependence of the parameter a, i.e. 
a = a' (T — T* fF ) , then the critical field has also a T— dependence, of the form: 

Pc= Pc o\T-T c MF \ 3/4 (24) 

where we have used the standard notations of the Landau theory of phase tran- 
sitions (see for instance [T3]). 

4 The electrical polarization of the anharmonic 
chain 

The polarization is proportional to the average displacement of atoms along the 
chain; according to (12), one obtains: 

(Q) = tanhwp + j — , p < p c (25) 

cosh top 

For small fields, 

(Q) =up(l + ••■) (26) 

and, comparing with the magnetization of the Ising model for small magnetic 
fields QI], 

M = l3He 2fiEl + 0((3 2 H 2 ) 
we can identify the coefficient u) as 

w ~ Pe 2f3c < (27) 

Qualitatively, (25) has a correct behaviour, in the sense that it is a monotonically 
increasing function of p, which saturates asymptotically. The same behaviour 
occurs in the ID Ising model in an external field. Also, the comparison with 
the Ising model recovers the temperature dependence of the numerical constant 
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w, entering in the solutions obtained in [S]. This recovery is essential in order 
to obtain the statistical mechanics of the systems under scrutiny. 

With (25), (27), we find for the susceptibility of the chain the following 
formula: 



2uj f uptaxihuip 



P<Pc (28) 



cosh 2 wp V coshwp 

5 The spatial array of Ginzburg-Landau chains 

Let us consider a planar array of N chains. The n— th chain (1 < n < N) is 
described by a Hamiltonian where is obtained from H c i, eq.(3), 

replacing u (x) — > u n (x) , p (x) — > p n (x) . Cyclic boundary conditions are im- 
posed, so that the chains n = 1 and n = N + l coincide. The interaction between 
chains n, n! is given by: 

H int = T ^2 D nn> [«n (x) - U n , (x)f (29) 

In a mean field approach, the Hamiltonian (29) is replaced by: 

-2qiD(u)J2 u n(x), * = 2,3 (30) 

n 

where g, is the number of near neighbors, in i dimensions; we shall focus here, 
however, only on the 2D case. The evaluation of the statistical physics of 
the 2D system reduces, throug the matrix transfer method, to the study of 
a Schrodinger equation describing a quartic oscillator in an external field T (u) 
[3]: 



' - 1 ~\A\u* + \bu*-T{u)u 



tpn (u) = e„Vn ( u ) (31) 



2m* du 2 2 1 1 4 
Measuring the energy in units \Vq\ — A 2 /AB and making the changes: 



(32) 



we get the following eigenvalue problem 
1 d 2 



-U-—. t — 4w +4u — 7 (u) u 
A du z 



■>jj n (u) = € n lp n (u) 



(33) 



where: 



0-1/2 1 (k B T\(\A\V /2 



(34) 
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In the classical paper of Bishop and Krumhansl [BK], the eigenvalue problem 
is solved numerically. Our contribution, to this point of the problem, is the 
following: using the analytical solution (anzatz) of Van der Straeten and Naudts 
[5], we shall obtain analytical formulae for the critical temperature of the 2D 
transition, and the critical behaviour of the planar array of chains. 

Let us now outline our approach to the 2- and 3D problems and describe the 
main results obtained. 

The mean field equation 

Replacing in (11) 



P = l(u) (35) 
we get the self-consistent equation: 

bx — x + tanha; — xtanh 2 x, b= — , x = uij (u) (36) 

For b < 2, the line described by the left hand side of the equation has a non- 
zero intersection with the curve described by tht right hand side, consequently 
a value (u) ^ does exist; so, an ordered state appears, meaning that the 2D 
system displays a phase transition. Consequently, 6 = 2 reprezints the critical 
condition, giving the expression of the transition temperature. 



Evaluation of the critical temperature T c 

According to eq. (27), 



k B T \k B T 
and the critical condition can be written as: 



ex p(i-^ cu o) ( 37 ) 



k B T 1 \k B T 

Putting 

2 



oj = — exp [ - cuq I (38) 



, T c< = Z (39) 
it takes the form of the transcendental Lambert-Euler equation: 

& = ^ (40) 

7 

Its solution is given by the Lambert function W [9] : 

f = W (f) ( 41 ) 
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The Lambert function has the following behaviour, for small, respectively large 
value of its argument: 

(42) 

So, for small values of the inter-chain coupling, proportional to 7, 

T c ~ pr^ (43) 
I In 71 

This result coincides to that obtained by Scalapino, Imry and Pincus [10] . using 
a different approach. For moderate values of the coupling, 

T c ~ 7 (44) 

The formulae (43), (44) confirm, qualitatively, the behaviour of T c , obtaind nu- 
merically in [15] . Actually, it corrects their interpretation for small 7 (considered 
by these authors to be exponential), and provides an analytical expression for 
T c , valid for any value of the inter-chain coupling. 



The critical behaviour of the order parameter 

In order to obtain the critical behaviour of the order parameter, we have to 
examine the self-consistent condition (36) for temperatures close to T c , where, 
putting lo{T c ) = w c , 

LU = lo c + Acj, Acj -C uj c (45) 

So, 

— = 2 - 2 7 Aw (46) 
CJ7 

Expanding the hyperbolic functions near the origin, we get: 

(u)~\T~T c \ 1/2 (47) 

showing a mean field critical behaviour, as we can expect, taking into ac- 
count the methodology itself. However, our approach is interesting, due to 
the fact that it also allows a simple, analytical determination of the coefficient 
of \T-T C \ 1/2 ■ 

The expression of the electrical susceptibility 

In the presence of an external electric field, the Hamitonian obtained through 
the transfer matrix approach is: 

hW> = --Li* + \a* + 1b„- - n«, n = e +7 <u> (48) 
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So, the ground state energy becomes: 

e (n) = e (0) - |a|IItanhwII (49) 
and the self-consistency condition takes the form: 

(it) = tanh (£ + 7 (u)) + u (£ + 7 (u)) [l - tanh 2 u (£ + 7 (u))] (50) 
which can be also written as: 



3£ 

(b- l)y-e = tanh y-y tanh 2 y, y = uj (£ + 7 (u» , e = -— (51) 

On this equation, one can easily see that the line defined by the left hand side 
has always an intersection with the curve defined by the right side, in a point 
having the abscise x > 0. Consequently, in the presence of an external field, 
(u) ^ 0, as expected. For small values of y, the self-consistency condition gives 
the following equation: 



-y 3 + (b - 2) y - e = 



(52) 



The physically interesting regime is still that corresponding to b < 2, so, putting 

AT 



= T 



we get the solution: 



(u) 



wf I 8 



1/3 



1/2' 



1/3 



11 1 



9e 2 ' 



1/2- 



1/3" 



-e 



(53) 

The exponent of the critical isotherm can be obtained immediatly, putting r = 
in the previous equation: 



(54) 



So, the exponent of the critical isotherm is S — 3, as expected from the general 
theory [Stanley]. 

However, the critical behaviour of the susceptibility can be obtained simpler, 
through the differentiation of (52). Keeping only the smallest terms, we have: 



dy 
de 



1 



(55) 



6-2 AT 

As this report is proportional to the electrical susceptibility, we find for its 
critical behaviour a Curie-type law: 



11 



The coefficient of the (T - T c ) factor can be easily obtained, using standard 
methods (see for instance [16], Ch.6). 



6 Conclusions 

In this paper, we have applied the transfer matrix method in order to study 
the statistical mechanics of ID and 2D Ginzburg-Landau systems - chains of 
anharmonic oscillators, or planar arrays of such chains, which can simultate 
ultra-thin films of ferroelectric systems. The starting point of our approach is 
a partially numeric / partially analytic expression for the ground state energy 
of a quartic oscillator in an electric field, eo (p), recently obtained by Van der 
Straeten and Naudts. In fact, the formula for eo (p) is mainly an ansatz than 
a rigorous result. Comparing the predictions of this formula, via the transfer 
matrix method, for the polarization of the Ginzburg-Landau chain in an electric 
field, with the small-field limit of the exact formula of the magnetization of the 
ID Ising model in a magnetic field, the temperature dependence of the ground 
state energy of the quartic oscillator in an external field is restored. Due to 
this fact, the expression of the ground state energy eo (p) can be used, through 
a mean field approximation, in order to find the statistical mechanics of the 
2D Ginzburg-Landau system. The method gives an analytical formula of the 
critical temperature of the 2D order/disorder transition, expressed as a Lambert 
function of the inverse inter-chain coupling. 

Through this approach, we are able not only to find the statistical mechanics 
of ID and 2D Ginzburg-Landau systems, but also to confirm the adequacy of 
the ansatz made by Van der Straeten and Naudts while proposing a formula for 
the ground state energy of a quartic oscillator in external field. 

Acknowledgement 1 The author thanks to CNCSIS and UEFISCSU for fund- 
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